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The microscopic approach was developed for obtaining of the free energy of a superconductor 
with help direct calculation of the vacuum amplitude. The functional of free energy of the spatially 
inhomogeneous superconductor in a magnetic field was obtained with help the developed approach. 
The obtained functional is generalization of Ginzburg-Landau functionals for any temperature, for 
arbitrary spatial variations of the order parameter and for the nonlocality of the order parameter 
and the magnetic response. Moreover the nonlocality of the magnetic response is the consequence 
of the order parameter's nonlocality. The extremals of this functional are considered in the explicit 
form in the low-temperature limit and in the high-temperature limit at the condition of slowness of 
spatial variations of the order parameter. 
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I. INTRODUCTION 



One of the important applications of condensed matter physics is description of thermodynamics and 
electrodynamics of superconductors. Several superconducting materials are spatially inhomogeneous anisotropic 
structures. In particular, heterostructures and superlattice, combinations of superconductors with the materials with 
different types of a conductivity. The superconductive and transport properties of the inhomogeneous systems is 
determined by the proximity effect, which results to modulation of the order parameter, and by Josephson effect 
P, Ql- At the present moment many various heterostructures have been made on the basic of the HTS. However 
the HTS can be considered as the heterostructures [s], where the conducting layers Cu02 play a role quantum 
wells, and the layers of atoms between them (the reservoirs of charge carrier) play a role quantum barriers. The 
value of the critical temperature can depend on the spatial structure of the superconducting system essentially. 
So, in the papers the model of the cluster crystal (the lattice is formed by metallic clusters) with the giant 

intensification of a superconductive pairing (hypothetically Tq > 150K) was proposed. Such crystal can be considered 
as a superconducting superlattice. The intensification of Cooper pairing can occur in a neighborhood of the special 
defects in a solid @. 

The modification of Ginzburg-Landau equation is necessary for descriptions of spatially inhomogeneous 
superconducting systems. So, the modified functional of free anergy is proposed for the description of the influence 
of the defects [9], where the term —'-fS{x)\A{x)\'^ is introduced. This term considers the intensification of Cooper 
pairing around the defect, moreover the area of intensification is much less than the correlation length Iq of the 
superconductor. For description of layered compounds of type HTS and MgB2 [3| the modified Eliashberg equations 
with consideration of anisotropy and Lawrence-Doniach functional ^] for quasi two dimensional superconductors are 
applied. The nonlocality causes additional complication in the description of the spatially inhomogeneous systems. 
The superconducting systems with a deep modulation of the order parameter must be considered as the systems of 
Josephson contacts. For the structures where d ^ A (d is the thickness of a contact, A is the magnetic penetration 
depth) the nonlocality is the determining factor of electrodynamics effects |lC)| - [l2| . However the effects of the 
nonlocality can be essential for the cases d <C A even (Tsl-fTsj . 

For the complete description of a superconductor it is necessary to know the functional of free energy i7(/3, A,A), 
where /? = y is the inverse temperature, A is the gap function (order parameter), A is the potential of microscopic 
magnetic field. Then the equations ^ = = are the equations of state of a superconductor. The basic method 

for obtaining of the equations of state is solution of Gor'kov equations or Bogoliubov equations with the equation 
of self-consistency for the order parameter [l6| A(r) and A+(r). However these set of equations is nonlinear. If the 
temperature is close to the critical temperature T — >■ Tc => A — > 0, then these equations can be represented in the 
form of series in degrees of A. Moreover, magnetic penetration depth A must be bigger then Pippard coherent length 
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Iq, hence the potential A changes a httle on a coherent length. As a result we have the well-known Ginzburg-Landau 
equations. The equations are correct for description of thermodynamics and electrodynamics of a superconductor at 
the following restrictions: 

1. The energy gap is much less than critical temperature. Then the parameter A{r, T)/Tc <C 1 can be expansion 
parameter. This means, that the equations are correct in the range T — > Tc or in the range H — > Hc2 (intensity 
of magnetic field is near the second critical magnetic field Hc2)- 

2. A(r,T) changes slowly on the coherent length ^T), which is size of a Cooper pair. 

3. Magnetic field H(r) = rotA(r) changes slowly on the coherent length, that is the magnetic penetration depth 
is A(T) ^ 2(0). This means, that electrodynamics of a superconductor is local. 

In the papers (itI . [18| the equations has been proposed, where the first restriction is absent. These equations 
were obtained from Gor'kov equations and they are the generalization of Ginzburg-Landau equations for the case of 
arbitrary value of A(r, T)/T. However spatial inhomogeneities are slow, order parameter and magnetic response are 
local. In the paper on basic Bogoliubov equation it was shown that if the gap A(r) varies with position, that the 
value of the gap in a point is determined by a distribution of gap in some neighborhood: A(r) = J (ir'Q(r, r')A'(r). 
This means that the order parameter is nonlocal. 

The problem of description of the superconducting phase can be solved by another method. This method is the 
direct calculation of the free energy. For it we must calculate the partition function Z 

Z = Sp (^ZoPoUW)) = ZoRiP), (1) 

where Zq is the partition function for a system of noninteracting particles, U{l3) = e+^'^^^^^^^e^'^^^^''^^ is the 
evolution operator in the interaction representation (it describes evolution of the system in imaginary time it — > r, 
T is the ordering operator in time), H = Hq + Hj is Hamiltonian of a system, N is the particle operator, fi is the 
chemical potential (replacement of Hamiltonian H in canonical ensemble by Hamiltonian H — fiN in grand canonical 
ensemble leads to the shift of reference of particle's energy from zero to Fermi surface: eikp) = 0), 

R{p) = {Umo = Sp (poUW)) - J2 —r / ■ • • / d^nSp (pof {i/i(T/) . . . i//(T„)}) (2) 

is the vacuum amplitude of the system, where Hi{(3) — e+''^^°^^^-'i?/e^'^^^''^''^-' is the interaction operator of 
particles in interaction representation. The averaging ()o = Sp(j)Q . . .) is done over ensemble of noninteracting particles. 

The partition function Zq can be found exactly for any system. If particles interact, then the situation becomes 
complicated essentially. Solution of the problem reduces to the calculation of a transition amplitude "vacuum- 
vacuum" ([21). The vacuum amplitude determine the internal energy of a system: 

t/ = -|ln^o-|lni?(/3), (3) 



and the grand thermodynamics potential: 



n = -ilnZo-^lni?(/3). (4) 



The potential ft plays a role of Helmholtz free energy with a reference of a particle's energy from Fermi surface. 
Therefore we shall call the grand thermodynamics potential by free energy for a brevity. 

The formulas pi4p is valid if the symmetries of ground state of the system with interaction |^'o) and without 
interaction |$o) are identical: ($o|^o) 0. This means that: l)potential of interaction is being switched slowly in the 
system in the ground state without interaction, 2) the ground state of the system with interaction is being obtained by 
continuous way from the ground state without interaction while the switching of the interaction {adiabatic hypothesis) 
(20I [2H. If the wave functions are orthogonal ($o|^'o) = - the adiabatic hypothesis is not valid, and the symmetries 
of the system with interaction and without one is different. This means, that the initial system without interaction 
suffers a phase transition stipulated by the interaction. For the system with the broken symmetry we can calculate 
the vacuum amplitude on the free propagators Go{'k,T) and use the formulas (|3I4I) . However we shall find the wrong 
thermodynamics potentials. This means that the state exists with more low energy than the found value. 

The method for direct calculation of the vacuum amplitude R{f3) has been proposed in the paper [2^. The concept 
lies in the fact that we consider electrons in a normal metal propagating in random "field "of fluctuations of the order 
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parameter Aq, where q is small wave-vector. The operator of the interaction of electrons with the fluctuations can 
be written as: 



E [^q^P,^-P_ + ^h^-P-Cp^] ' (5) 
p 

where p_|_ = p ± q/2. The correction to the thermodynamics potential is represented via the vacuum amplitude R 
using Wick theorem: 

1 /X" 1 /X" 

An = -TlnR{l3) « -T[R{I3) - I] ^ -L f [ dT2{fr{H^nt{T^)H^nt{T2))) +■■■ 

= -T I dn I dT2|Aq|2 VGo(p+,Ti-T2)Go(-P_,Ti-T2) + .... (6) 
Jo Jo T-> 



The correction Ail is represented via the free propagators Go of normal state only - we consider normal metal at 
T > where the fluctuation sourse of Cooper pair (j36p acts. As a result we have Landau expansion: 



q 



a(T)|Aqp + ^|Aq|4 + 7g'|Aq| 



(7) 



where a{T) cx {T — T),b,^ are expansion coefficients. As for calculation of the correction Ail the free propagators Gq 
of normal phase is used only, the condensed phase is considered as fluctuations against the background of the normal 
phase. This means, that we can obtain the limit f2(T — > T) only. Moreover, the hamiltonian ^ doesn't contain any 
parameters of matter (for example, interaction constant between fermions). 

In the presented paper we have developed the microscopic approach for a finding of the functional of free energy 
of a superconductor J1(A,T) with help of direct calculation of the vacuum amplitude. This value is calculated 
on the dressed one-particle propagators. The propagators is dressed due to the interaction of free fermions with 
the fluctuations of pairing by the method presented in [23, 24]. In the sections IIIII and IIVI using the developed 
approach of microscopic description of superconducting phase we obtained the functional of free energy of a spatial 
inhomogeneous superconductor in a magnetic field Jl(Aq, a(q), T). This functional generalizes Ginzburg-Landau 
functional for arbitrary temperatures, arbitrary spatial variations of the order parameter and for nonlocality of order 
parameter and magnetic response. 

II. THE METHOD OF UNCOUPLING OF CORRELATIONS FOR CALCULATION OF THE VACUUM 

AMPLITUDE. 

A. Normal and anomalous propagators. 

Let we have the system from N noninteracting fermions in volume V at temperature T. For description of the system 
we must use Matsubara propagators, where time t is complex: t —ir, r e [0, (3]. In ideal Fermy gas propagation of 
a particle with momentum k, with energy e « ^^iJ-dkl — kp) counted off Fermy surface (we are using system of units, 
where h = ks = 1) and wiyh spin a is described by the free propagator: 




r > 
r < 



*Sp{poCk,,(r2)C+^(ri)} 

1 n=+oo „^ 

G(k,T) = - ^ G(k,a;„)e-^'^"^ G{k,Lo^) = - J G(k, T)e^""-dr, 



1, T > 

0, T < 



(2n+ l)7r 



where 



r< n \ * • ^'^^ + e . Aq Bo 
Go(k,w„) = ^ TTT^i-r. — 72 2^^- M+*^ Tm' 
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Cjj ^(r) and (t) are operators of creation and annihilation in Heisenberg representation, po is density matrix of 
noninteracting particles: 

Po = exp <^ — - — ^ = exp <^ S . (10) 

Now let an attracting force acts between particles. Hamiltonian of the system is 



^0 + i?/ = E E ^ (k) c^A'^ + ^ E E vi-ik^-k9^''^\SA^~k,, (11) 

a k "■'T k,l 

The force is described by matrix element of interaction potential: 



(1, -\\V\k, -k) = Xwiwk < 0, = { J; 'ffcj 5 ^ (12) 

1, a^j 



moreover interacting particles have opposite spins: ga-y ~ ^ q a — j 

The interaction (jl2p leads to Cooper instability. Let the additional particle with momentum k, lu propagates through 
the system of identical fermions. Some pair of fermions decays in components with momentums — k, —uj and k, w with 
amplitude iA+. Second particle of the decayed pair is in a state of the additional particle (k, lo) and it is identical to 
the additional particle. The second particle propagates through the system further. First particle of decayed pair forms 
bound state with the initial additional particle with amplitude — lA. Anew formed pair replenishes the condensate of 
pairs in the system. Thus, the dressed propagator Gs considers the interaction of a particle with the fluctuations of 
pairing. Intensity of the interaction is the amplitudes —iA and iA~^ . Then we can write the mass operator for such 
process (FigH]) as [H, ill 

-AA+ 

- S(k,a;„) = (-A)zGo+(-k,^„)(-A+) = . (13) 

iLOn + e(k) 




Phc. 1: The diagram for the mass operator E describing interaction of a fermion with fluctuations of pairing. 

Then it follows from Dyson equation, that a dressed propagator has the form: 

11 i 
iE Gs(k,w„) 



Go Gs ' icjn - £(fc) - E(k, w„) 

_ iUJn + £ _ As Bs 

~\iWn? - E'^{k) ~\uJn- Eik) \uJn+E{ky 



(14) 



where E ~ y^e^{k) + AA+. It can be written with help of a total definition of Green function in (k, t)-space: 

G.(k,.)^.-n)J-"^>^^>^)^^^^^ 

\ iSp[goC+Jn)G^jT,)}, r<0 

= -iOriajAse-''^ + gsBse""^) + lO-Ags Ase'""^ + gtBse""^), (15) 

where 

^s-\{^ + ^)^ Bs-\{l-^), 9t-^^^. ^s-^,- (16) 
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^0 is the density matrix of noninteracting quasi-particles: 



00 = exp ■ 



T 



(17) 



Let's introduce the designations: — GS = AF^, — = A'^F. Then Dyson equation can be rewritten in a form 
of Gor'kov equations: 



{iuj„ - e)G + AF+ = i 
(iion + e)F+ + GA = 0. 

The expressions for anomalous propagators follow from Gor'kov equations: 

-iA+ 



F+(k,w„ 



F(k,c^„) = (^+(k,^„))+ = 



iA 



{i^r^Y ~ E^{ky 



(18) 
(19) 



(20) 



We can write the anomalous propagators in (k,t)-space and in a form of a vacuum average of creation and annihilation 
operators: 



i^a7(k, r) 



A+ r *Sp{?oG+j^^(r2)C+ Jn)}, r > I 

VA+A \ *Sp{?oG+ Jri)G+j^^^^(r2)}, r < J 

V A+A 

A r -*Sp{^oGk^(r2)G_k^(n)}, r>0 



TA+A I -*Sp{?oG_k,^(Ti)Ck^^(r2)}, r < 
A 



(21) 



(22) 



we can see, that the existence of nonzero anomalous propagators F and F'^ means breakdown of global gauge symmetry 
in a system, that is number of particles is not conserved in the course of existence of a pair condensate. 



B. Kinetic energy and entropy. 



In order to calculate a free energy it is necessary to know kinetic energy of particles of a system, energy of their 
interaction and entropy of collective excitations. Average kinetic energy of all particles of a system is 



{W) = -2z^G(k,r->0-)e(fc) = -^ lim_ XI G(k, w„)e-'""*e(fc) 
k k «= 00 

/OO 
{g-A + g+B)ede, 



n— +00 



(23) 



where — ^^'^ density of states on Fermy surface. Since the interaction (|T2l) exists only in the layer 

—ujjj < e{k) < iojj, then we can suppose that 



g 



9o ; > 
9s'-: < 



5+ = 



g+: \eik)\ > ujd 
g+; \e{k)\ < ud 



Hence, one may write expression for kinetic energy separating normal and superconductive parts as 

r° f 3\e\\ f^" / 0E\ 

(W) = W„ + Viyp J tanh - Vvp J tanh J -^de. 



(24) 



(25) 



If we suppose that A = 0, then we shall have W — Wn- 
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At temperature T 7^ a gas of collective excitation exists - boholons with the spectrum E = y^s^{k) + A^. Since 
boholons are product of decay of Cooper pairs on fermions, hence occupation numbers of states by boholons are 

/sW = ^4(l-tanh(f)). (26) 

Then entropy of a system is 

S = -2^[/(fc)ln/(fc) + (l-/(fc))ln(l-/(fc))] 
k 

[Js In /s + (1 - Is) ln(l - fs)] ds + 2Vvf / [/o In /o + (1 - /o) ln(l - /o)] de. (27) 

Here we separated the normal part again, where /o = (e^l*^! + 1)"\ so that 5 = 5„ at A = 0. The muhiplier 
"2"appeared as result of summation over spin states. 



C. Vacuum amplitude. 

In the previous subsections we considered the interaction of particles with fluctuations of pairing, and we found, 
that the state of a system described by the density matrix ^0 has another symmetry in comparison with the initial 
state po- The effective Hamiltonian of the system of quasiparticle is 

+ ^/ - E E ^ (k) C^Ac + ^ E E n,-Lk,-kff«7C+i_^C+„Ck^„C_k,, (28) 

a k " T k,l 



The vacuum amplitude of a system can be written in the form: 



R{I3) = ([/(/3))o - Sp [goUiP)) = ^ / rfn . . . / dr^^Sp (goT [h,{ti) . . . il/(T„)}) , 



(29) 



where Hj{t) = e'^'^^" Hje^'^^" is the interaction operator of particles in interaction representation. The averaging 
Sp ^goC^(/3)^ is done over ensemble of noninteracting quasi-particles. It can write the expended expression for the 
vacuum amplitude: 

= 1 + W r^^iEE (-^^l,-l,k,-k) Sp{?oC+ ^(n)C+Jn)Ck,,(n)C_k,,(n)} 

■ •'0 a,f k,l ^ ^ 

""J ""J a,7 k,l ^ ^ "'^7'k',r ^ 

xSp{eoC+,^/r2)C+^X-2)Ck^„,(r2)C_k^,,(r2)C+^^(rOq+JrOCk^Jn 

+..., (30) 

where we took into account that Sp{^o} = 1- In order to calculate ([50]) approximately we can uncouple correlations 
by the following way taking into account anticommutation of operators C and C"*": 



m « 1 + ^(-1)'^ / '^^i E E (-^^l.-l,k.-k) Sp {?oq+Jn)C+ ^(n)} Sp {doC_],Jr,)C^Jn)} 

a,7 k,I 

"'O ""^ a,7 k,l ^ ^ "'.7'k',l' ^ ^ 

xSp{aC+^,(r2)C+,^,(r2)}sp{^oC_k^,,(r2)Ck^,„,(r2)}sp{^oC+Jn 

+ ... = l + Ei + ^i?? + ...=exp(i?i) (31) 
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Let's take into account that our approximation is analogous to Fock approximation for normal processes. A decay 
of quasi-particles is absent in Hartree-Fock approximation, hence the amplitude of pairing is real A = A+ in the 
momentum space. Hence F = — F+. Then R{t) can be written as 



■■n k.l 



T7 E f J^L-Lk.-k) n\ r ^ 0-)F(k, r ^ 0-)/3 



lni?(/3) =i?i(/3) 



V 

kj 

^ n— +00 ^ n— +00 

n— — 00 n— — 00 

= /"J^ tanh Arf, tanh Arf,. (32) 

The procedure of uncoupling of correlations and the contribution of the anomalous process in the vacuum amplitude 
can represented graphically (Figl2]). In order to calculate a contribution of interaction in free energy we can use the 
formula 

^ A 1 \ n— + CXD 

n— — oo ]^ n— — oo 

= XV.I tanh (^^) j'^J tanh Arf^. (33) 

If we suppose that A = 0, then we shall have Q,x = 0. 




Phc. 2: The uncoupling of correlations in a vacuum amplitude for the process of scattering of two fermions by the potential of 
interaction Vj _j _^ from the initial state k, a and — k, 7 to the final states 1, a and —1, 7. As a result we have the anomalous 
transition amplitude "vacuum- vacuum". 



D. EYee energy. 

Starting from the above found results we can write the expression for free energy of a system: 

„ . n=+oo 

^={W)--S + nx^~^ lim V V G{li,uJn)e-'^-'e{k) 

n— — 00 

^ , ^ n— +00 ^ 71— +00 

+ -^[/(fc)ln/(fc) + (l-/(fc))ln(I-/(fc))]--5]z«,- ^ F(k,c.„)E^'^-« E (34) 

k n— — 00 ]^ n— — 00 
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The observed value of A must minimize the free energy: 

— = 0=^H)A=— ^ ^ w,F{k,u;,,). (35) 

The functional (j34p can be written in quadratures: 



n— — oo 



de 



+ ^'^F / " [/5 In /s + (1 - fs) Hi - /s) - /o In /o - (1 - /o) ln(l - /o)] de 
+VVF9 I " tanh ^de f " tanh ( ^) ^ds, (36) 



^2 J 2E J_^^ \2 J2E 

where f2„ is the energy of a normal phase, g = Xvp is the effective interaction constant, V is volume of a system. If 
to suppose A = 0, then we shall have = ri„. The equilibrium value of A is determined by balance of kinetic energy, 
the entropy term and energy of interaction, g can be expressed via critical temperature Pc with help of the equation 
^{Pc) = as following: 

Let's consider a low-temperature limit of the free energy (^5)1 : A/3 ^ 1 at T — 0. This means, that the value A — Aq 
can be parameter of expansion, where Aq — A(T = 0) is equilibrium value of the gap (amplitude of pairing) at zero 
temperature. Due a rapid convergence of integration elements in p6p . the limits of integration can be — cx), +oo. Then 
a low-temperature expansion has a form: 

= a. + "t^(ao(T)-K6o(r)A + doA2), (38) 
where coefficients of the expansion for the approximation g <C 1 A ^ are 



boiT) = -2i/i.Ao + i/F^STrAoTe"^, do = (39) 

The equilibrium value of the energy gap is A(r) — — Aq ^1 — ^ 

Let's consider a high-temperature limit of the free energy: A/3c <C 1 at T ^ Tc- Expansion in powers of A gives: 

n^n„ + V (a{T)A^ + IbA^ + IdA") , (40) 



2 3 



where the coefficients of the expansion are 



m ^-^c , 7C(3) 

aiT)^.F^-, b = .F^ (41) 

/52.31C(5) ^ „\ 1 
d = vf{ + 4.83 ' 



4!T4 ■ 

This expansion has a form of Landau expansion of free energy in powers of order parameter. 



III. FREE ENERGY OF THE SPATIALLY INHOMOGENEOUS SUPERCONDUCTOR. 

In the previous sections we supposed, that amplitudes of pairing A and A+ don't depend on spatial coordinates. 
This takes place in interminable, homogeneous, isotropic and isolated from external fields superconductor. However 
in a total case these conditions are not realized. For example, in a sufficiently strong magnetic field the inclusions 
of normal phase can exist in volume of a superconductor. Another example: the contact of a superconductor and a 



9 



normal metal. In this case the order parameter is suppressed in the boundary layer of a superconductor, however 
it appears in the boundary layer of a normal metal. We understand a spatially inhomogeneous superconductor as a 
superconductor with spatially inhomogeneous distribution of the order parameter: A(r) and A+(r). 

In the previous sections we considered the pairing of fermions with the opposite momentums: k and — k. Now let 
fermions can pair with arbitrary momentums k + q and — k, where q passes through all the vector space. The potential 
of particles' interaction in the given states is Vj_|^q _j i5;_|_q \^ = ^WkWi, A < 0. For convenience let's pass to the 
reference system of center of mass of a pair, where the momentums of particles are equal by modulus and are directed 
oppositely: 

_ TO2P1 - TO1P2 _ , ^ q 

Poi — ^ + o' 

mim2 2 
rn^p^-rn^^_^_^ (42) 
mim2 2 

because mi — 7712 = m. Hence the hamiltonian of the system has the form: 

Se«-Ei:i:Ki'+f)c;,,Ai,„ 

q " k 

+ ^EEE^«.^lq,S"a„^k+a,„^_k_av (43) 

q k,l 2.2' 

In the system with such hamiltonian the Cooper instability to takes place either as for the system with the hamiltonian 
However now the amplitude of pairing is function of q: A — > Aq. If the energy gap is the function of a wave 
vector A = A(q) then the energy gap is the function of a radius- vector: A — A(r) because the Fourier transformation 
realizes the one to one correspondence between functions defined in the q-space and in the r-space. Taking into account 
the reality of the order parameter in the momentum space Aq = Aq and requiring the identical dimension of the 
order parameter in the q-space and in the r-space we can write the Fourier-transformations as: 

A(r) = ^ A(q)e'q'- ^ J Ai<i)e^^^ d'q, 

A+(r) = ^ A(g)e-qr ^ (2^ / ^(q)^''"""'^'?' (44) 
A(q) ^ IJ A(r)e-*qi'd3r = i J A+(r)e'q^d3 




Phc. 3: The diagram for the mass operator E describing an interaction of a fermion with fluctuations of pairing in the spatially 
inhomogeneous system. 

The mass operator for the process of interaction of a particle with the fluctuations of pairing is shown in Fig|31 In 
analytical representation it has a view: 

—A A"*" 

- Sq(k,c.„) - (-Aq)zG+ (-k - |,c.„) (-A+) = (45) 



iujn -I- e ( k + ^ 



where the free propagator Go is 



n 1 . iuJn + eq{k) 

Go = 7 r = *7 TT7T- (40) 
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Then from Dyson equation we can obtain the dressed propagator: 



1 1 .„ „ . iw„ + eq 

Go Gs yi^^nY - 

where i?q is the specter of quasi-particles in a nonhomogeneous system: 



El^ = e\ + |Aq| 



£q = £ I k 



e(fc) 



2m 



|q| < fcF- 



(47) 



(48) 



Dyson equation can be represented in a form of set of Gor'kov equations, from where the expressions for anomalous 
propagators follow: 



(zw„ - £q)G + AqF+ = i 
{iujn + eq)i^+ + GAq = 



^^q(k,L^„) = (F+(k,c^„) 



(49) 



If to suppose (7 = 0, then we shall have the expressions (I18m9p . 

Now let's suppose that A(r) changes very slowly on a coherence length 1{T) which characterizes a size of Cooper 

pair. Then we can suppose £ ^k + = £(k) in the specter of quasi-particles, such that i?q ^^e'^{k) + |Aqp. 
However we must keep £ (k + ^ ) in numerator of the expressions p6)) and (|47)) for G. Hence the normal propagator 



5+A(k + q)e-^q^ + 5qi?(k + q)e^q^ 



has a form: 
Gq(k,T) 
where 

^(k+l) «Aq(fc) 
The anomalous propagators are 



.gqA(k + q)e-^q^+5+B(k + q)e^qn , (50) 



1 kq 



B(k+|) «i?q(fc) 



1 kq 

2Eq2rn' 



(51) 



Fa^(k, q,T) 



Aq 

2Ea 



(52) 



We can see, that in the approximation of slowness of changes of A(r) the anomalous propagators depend on q by 
means of A(q) only. 

Kinetic energy of a system is determined by the following way: 



(»') 



0" 



q k 



+ Vvp 



tanh 



tanh 



(PEc^\ 1 



m 
2 



tanh 



2 



Er, 



de 



de. 



(53) 



We can see, that the term, which is proportional to g^, is added to the kinetic energy (with the replacement 
E — > i?q = ^Je^ik) + |A(q)p) . In the approximation of slowness of changes the expressions for entropy and vacuum 
amplitude coincide with the expressions (P71) and accordingly, however it should be written A(q) instead of 
A — const. Then we can write the free energy: 



r! = a. + 5] f^(Aq) + Vvp'^q^ J2 n 
CI n •'-'^D 



tanh 



tanh 



V 2 j i?q. 



(54) 
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where ri(Aq) coincides with the expression where the replacement A — )■ A(q) was done. 
Expanding the free energy (j54p in powers of A we can obtain the expression: 



(55) 



where the coefficient 7 is 



1 = 



487r2T2 



(56) 



where Iq is a coherence length at T = (Pippard length). The expansion ([55]) has a form of Landau expansion of free 
energy in powers of order parameter at the condition qIq <^ 1. We can see, that a spatial inhomogeneity increases the 



free energy of a superconductor. Hence in most cases we can be restricted by the term - 
more fast changes of A increase the free energy essentially. 

Let's pass from momentum space to real space using the expressions (|44l) : 



q in the expansion, because 



AqAqd^g 



1 Aq 1 1 A(r)e-'q'-dV d\ ^ ^ J A(r) J 



A{q)e-'^''d^q 



d-'r 



1/2 



J A(r)A+(r)d3r,(57) 



qAqqAqd^g^ J qAq 
qAiqje-'^^^d^q 



e-^^T^(-i)^A(r)d\ 
or 



d^q 



V 



dr 



-A(r) 



d_ 
dr 



A+{r)d^ 



(58) 



For the term A'* and terms with more high powers the situation is more difficult. This is because a square of a Fourier 



transform is not equal to a Fourier transform of a square: (-^ J A{r)( 



r) ^ ^jA'^{r)e''^^d^r. Apparently 



this fact is the manifestation of a nonlocality of the order parameter in zero magnetic field described in the paper [19|, 
where the value of gap in a point is determined by a distribution of the gap in the some neighborhood of this point: 
A(r) = / dr'Q(r, r')A'(r), where Q{r,r') is the kernel in coordinate space. If we neglect of this correlation then we 
shall obtain the expansion of the free energy in powers of AA"*" in real space in the following form: 



— 



a(T)|A(r)|2 + -|A(r)r 



d_ 
dr 



A(r) 



d^j 



(59) 



This expansion coincides with Ginzburg-Landau expansion in zero magnetic field. As a result of minimization of the 
functional (|59p !£■ = with a consideration of the corresponding boundary condition we obtain the energy gap as the 
function of coordinates A = A(r). It prove our statement about that the pairing of fermions with the momentums 



and — k 



q 

2 



(or k + Y and — k — ^ in system of center of mass) results to the spatially inhomogeneous 



order parameter if q passes through all the vector space. The physical interpretation of this fact lies in the following. 
Let's assume that we have the homogeneous superconductor with spatially homogeneous distribution of the order 
parameter. This distribution corresponds to a minimum of the free energy. Then we shall create an nonhomogeneity 
in the distribution of the energy gap by some way. This nonhomogeneity increases the free energy according to (j54p or 
(jSOp . This means that we is doing some work on the system. This work is spent on infusion of an additional momentum 
to each partner of a pair: k k + ^ and — k — k — 3. 

Minimization of the functional (IsSi) gives (at the condition q <C l(T)): 



A2 

^q 



-^(l + lq^)=A^{T){l-P{T)q^) 



A2(r) 
i + ;2(T)(j2 



Q(q) 



1 



i + z2(r)g2/2' 



(60) 



where Q(q) is the kernel in momentum space, 1{T) = lo/\a{T)\ is the coherent length at temperature T. At the 
critical temperature T — >■ Tq we have l(Tc) — 00. This means that at high temperatures the nonlocality is caused by 
fluctuations in the critical region. 

Let's consider the low-temperature expansion: 



(MT) + &o(T)Aq + doAq + 79' A^) , 

q 



(61) 
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where coefficients aQ{T), bQ{T), do are determined by the formulas ([5^ . and coefficient 7 is determined by the formula 
psp . The observed configuration of the order parameter Aq minimized the free energy: 

g . „ ^ ,„(r) + 2*A<, + -„=A<, . ^ A<,(T) . . ^ Q . (62) 

Thus the energy gap (the amplitude of pairing) Aq is the nonlocal order parameter with the correlation length Iq. 

The dependence of the interaction constant on coordinates A(r) = J^q^i^)^^*^^ the direct cause of the spatial 
inhomogeneity. 

Heff - EEE^.(k + q)C^i+q.C^k+q.. 

q " k 

+ 2l7^5IE^(q)^fc^'^'"7C-l-q^^r+qa^k+q,a^-k-q,r (^3) 
Hence the free energy of a spatially inhomogeneous superconductor has the form: 

C[ n— — 00 

+ E E I [/q 1^ /q + (1 - /q) ^^^(1 /q) " /q /q " " /q) " /q,a) 

q k,(|e(fc)|<a;D,|k|~fcF) 



n— + CXD ^ n— + 00 



q k n— — 00 ]^ n— — 00 



the observed value of Aq.a must minimize the free energy 



The high-temperature expansion for the free energy can be obtained by directing the temperature T — >■ Tc(q), 
where the temperature Tc{ci) can be obtained from the equation (j68p assuming A = and £q = e: 

where g{ci) = A(q)i^F- Then Aq — > and the expansion of the expression (I64p in powers of Aq in the long-wave limit 
has the form: 



n V rc(q) ^ 2 87r2r2(q) q ^ 487r2T2 (q) q 



(67) 



Moreover the correlation length can be introduced formally: 7(q) = lypl-oi^)- the r-space the coefficients of the 
expansion is functions of coordinates a{T, r), 6(r), 7(r) and they can are obtained by the microscopic way only. 

The value Ao(q) — Aq(T) can be the parameter of the low temperature expansion, where the value Ao(q) can be 
obtained from the equation (|68l) assuming T — and eq = e: 



Then the low-temperature expansion has the form (j6ip in the long- wave limit, where the coefficients 
ao{T,q),bQ{q){T),dQ determined by the formulas ([M)) with the replacement Aq — ?> Ao(q), and coefficient 7 is 
determined by the formula ([55)) with the replacement Iq — ^ ^o(q)- In the r-space the coefficients is functions of 
coordinates ao{T,r),bo{r){T),'-f{r) and they can are obtained by the microscopic way only. 
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IV. FREE ENERGY OF THE SPATIALLY INHOMOGENEOUS SUPERCONDUCTOR IN A 

MAGNETIC FIELD. 



In this section we shall generalize the previous results for the case, when a superconductor is placed in a magnetic 
field H(r) = rotA(r). Our aim is to obtain the functional of free energy (A(r), ^A(r), A(r)) , which is correct for an 
arbitrary value of the relation A(r) /T, for an arbitrary scale of a change of A(r) in comparison with a coherent length 
Z(T), for an arbitrary value of a magnetic penetration depth A(T) in comparison with a coherent length Iq (nonlocal 
electromagnetic response). Thus, the all three restriction on Ginzburg-Landau functional, described in section U are 
excluded. 

Let the exact microscopic magnetic field with a potential A(r) and an intensity H(r) exists in the given point of 
the superconductor: 

A(r) = J2 aqe''^"" H(r) = rotA(r) =i^qx aqe*^''. (69) 

q q 



Moreover a_q = aq follows from the reality ofA. Energy of the magnetic field is 



Wf = ^l |H(r)|^ d'^-^T. (^'(aq • 4) - (q • ^q)(q • 4) 



(70) 



A magnetic field affects on a superconductor essentially. In the first place, a distribution of order parameter becomes 
inhomogeneous. As it was shown in the section lllli an inhomogeneity leads to the growth of momentum of each 
element of a pair: k — > k + ^ ^ ^ ^ moreover the order parameter depends on the momentum A = A(q). 
In the second place, an ordinary momentum must be replaced by a canonical momentum: k+^ — >k+§ — |aq 

and — k — — > — k ~ ^ + f^^-qi moreover the order parameter is function of the momentums: A = A(q, aq). For 
convenience let's pass to the reference system of center of mass of a pair, where the momentums of particles are equal 
by modulus, are directed oppositely and real: 

rn2Pi - miP2 , , q e . , +n 
Poi = k + - - — aq + a+ 



77117712 2 2c 

| + -,.q^«q, 



m2Pi - miP2 1 q , e , , 

Poi = = -k- + T^(aq + a+) (71) 



r71lTO2 

because mi — m2 = m. The the hamiltonian of the system has the form: 

Heff = EEE^?(k + q-7aq)c+^ ^^Ck+q_,aq,o 

q " k ^ 

+ ^ E E E ^(q)«^'=^'5a,C+ Ck+q_,aq.aC^-k-q+taq,v (^2) 

q " k,i 

if A(q) < 0. The hamiltonian (|72|) is the generalization of BCS hamiltonian to the cases of the dependence of the 
interaction constant of particles on coordinates and is the generalization of the hamiltonian ()6ip to the presence of a 
magnetic field. 

The system with the hamiltonian (|72p is unstable regarding in the Cooper pairing. The mass operator for a process 
of interaction of a fermion with a fluctuation of pairing is shown in FigU In analytical representation this diagram 
has the form: 

where the free propagator Go is 

I iw„ + £ (k+ § - ^(aq + ai) 
Go-- 7 (74) 



+ 2 - ^(aq + a^) iiLOnV - k + § - ^(aq + a+ 
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Phc. 4: The diagram for the mass operator S describing the interaction of a charged fermion with fluctuations of pairing in a 
spatially inhomogeneous system situated in magnetic field with a potential a(q). 



Then from Dyson equation we can obtain the dressed propagator: 



1 1 iw„ + e(k+§-^(aq + a+)) 

— = zSq=^G5-z \. (75) 

Go Gs ^ [lUJnr - Eq a 



where E is the specter of quasi-particles in the inhomogeneous system situated in magnetic field: 

2 2c 



= (k + ^ - f (aq + a+)) + | Aq.aP 



£(k+^--(aq + a+)j «£(A;) + ^ ^ = £q,a, |k| ~ A:^ , (76) 

where we have introduced the notations i?q,a and £q,a for convenience. Then Dyson equation can be represented in 
the form of Gor'kov equations. From these equations the expressions for anomalous propagators follow: 

{iujn - £q,a)G + Aq.aF+ = i ^q,a = (»:m^~e^ 

(77) 

{iujn + £q.a)F+ + GAna = ^ \+ ^^qa 

%a-(i^q,a) = (-.)^-B^.a 

If to suppose q — K a = 0, then we shall have the expressions p8ll20p . 
Free energy of a superconductor is the sum of the following terms: 

Q = nn + Ws-^Ss + nx + Wfia.), (78) 

where Ws is the kinetic energy of fermions of a system in superconductive phase, Ss is the entropy of boholons, 
= — ^ln^(/3) is the energy corresponding to an interaction, W/(a) is the energy of the magnetic field ([70)) . The 
expressions for Ws, S and are obtained from the expression (IM)) by the replacement A — >■ Aq.a, E i?q.a, 
e — > £ f k + ^ — ^(aq + aq) j , / — > /q,a- Hence the free energy of a superconductor 



IS 



^ = -5 EE E G(k+|~^(aq + a+,..„)e— s(k+|-^(aq + a+)) 

^ k fi'——oo 

+ E E I [/q.a In /q,a + (1 - /q.a) ln(l - /^a) - ^a In /^,a - (1 - /^,a) ln(l - /^,; 

q k,(|e(fc)|<a;D,|k|~fcF) 

-E-^E^fe« E %a(k,cj„)E^*^-« E %a(k,cj„) 

q k n— — 00 ]^ n— — cxD 



"Stt 



^(g2(aq.a+)-(q.aq)(q.a+)) (79) 

q 
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where 

fiaik) = ^pEci_aW + ^ ' ^q^al^) - ^^|eq,a(fc)l ^ ^ (^^^ 

are occupation numbers of states by boholons, moreover with help of the term /q ^ the normal part of the entropy is 
separated, such that ris(A = 0) = 0. 

Unlike Ginzhurg- Landau Junctional the obtained functional of free energy (j79p is correct for an arbitrary value of 
the relation A(T)/T, for any temperature, for arbitrary spatial variations of order parameter A(r) in comparison with 
a coherent length l(T), for an arbitrary value of a magnetic penetration depth A(T) in comparison with a coherent 
length Iq. The observed value of Aq.a must minimize the free energy: 

^ = =^H)Aq,a = ^5:"Xr-^^q.a(k,-J (81) 

If we suppose a = and neglect by the dependence on q then we obtain the expression p5l) . However the obtained 
expressions is complicated for analyze. Let's suppose that A(r) changes in space slowly. Then it is necessary to expand 
the expression ([79)) in degrees of ^ — ^(aq + aq) keeping terms which are proportional to the vector in second degree 
only. In order to simplify our consideration let's suppose that the parameter of interaction A(q) = const. This means 
that a spatially inhomogeneity of the order parameter is caused by the magnetic field only. Let's consider the high- 
temperature limit of the free energy: A/3c <C 1 at T — >■ T^. This means, that the expression ([7^ can be expended in 
degrees of Aq^a: 

fi = r!„ + y^(^a(r)A2j_a + ^6A^^a + 7(q-^(aq + a+))'A^, 

+ ^ E («'(^q • 4) - (q • ^q)(q • 4)) ' (82) 
q 

where the coefficients a{T), b are determined by the formula (|4T|) . and the coefficient 7 is determined by the formula 
([56| . The expansion (|82|) has the form of Ginzburg-Landau expansion of the free energy in degrees of the order 
parameter. Observed configuration of the order parameter Aq_a and the magnetic field a(q) minimizes the free 
energy: 



sn 



e , 



^ = ^ «(T)Aq^a + &A^ a + 7 (q - -(aq + a+)j Aq.a = (83) 



sn 

^ = ^ j+(q) - e7qA^,a " 7-A^.a(aq + a+) (84) 

— = ^ j(q) = eTqA^j^a " 7-A^.a(aq + a+) (85) 
where j(q) is Fourier component of a current: 

J(i") = £:^°*H(r) ^ j(q) = -£:q X q X aq = (q(q- aq) - aqg^) 



j+(q) = -^q X q X a+ - -— (q(q • a+) - a+g^j . (86) 

From the equations (|84l85p we can see that j(q) = j^(q). The from the equation (l86l) we can see that aq = aq. 

If the superconductor is simply connected (without holes, vortexes) we can pass to the transverse gauge of magnetic 
field q • aq — 0. Then the current is j(q) = ^a(q)q^. In transverse gauge the functional of free energy has the form: 



The equations of extremals are 



a(T)Aq,a + ?>A^ a + 7 (^q' + ^^q) ^^-a ^ (^8) 
j(q) = -87^A2jaaq. (89) 
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The closed currents q • j(q) = screen a magnetic field in a superconductor. The currents is analogy to molecular 
currents of Ampere. Their resulting results in observed magnetic effects. From (|89p one can see, that the value 
Q — — 87^Aqa^ is the kernel of magnetic response. The order parameter is function of q and a(q): 



A4.a(7^) = 1 - ri + — 4 H 7 / , , NN . (90) 



where ql{T) <^ 1. From the formula (^0)) one can see, that the kernel Q is a function of the magnetic field. Hence the 
electrodynamics of a superconductor is nonlinear. The correlation length 1{T) in the kernel of the magnetic response 
considers both the own nonlocality of a superconductor and the nonlocality caused by fluctuations in the critical 
region. If to suppose A = const at given temperature, then we shall obtain London equation: 

j(q) = -87-A2(T)aq ^ -^^a(q) ^ k^T) ^ f ^ ^ , (91) 

where A(T) is the magnetic penetration depth in a superconductor. 

For research of the nonlocal characteristics of the functional of free energy ([7^ let's consider a low-temperature 
limit A/3 ^ 1 at T — > 0. The value of a gap is close to the value at zero temperature A(T) < Aq. Moreover, magnetic 
field is weak, such that it changes a value of gap lightly, that is the magnetic field is much smaller than critical field 
H <C He- Either as above, we assume that a change of a gap in space is slow. Starting from aforesaid and using the 
expansion (j38p we obtain the free energy: 

17 = a. + V^J] (^ao(T) + 6o(r)Aq,a + doA^j a + 7 (q' + + ^ E '^'4' (^2) 

where coefficients af){T), bo(T), do are determined by the formulas p9p . and coefficient 7 is determined by the formula 
([56| . the magnetic field is considered in the transverse gauge q-aq = 0. The observed configurations of order parameter 
Aq.a and magnetic field a(q) minimized the free energy: 

1^ = ^ bo{T) + 2doAq^a + 7 (^q" + ^^q) ^^'^ = (93) 

AO p2 

-^=0 ^ j(q)=-87-A^.aa(q). (94) 
oa c ^ 

If in the equation (|94[) to assume A = const at given temperature, then we shall have London equation again: 



"0 

47rA2(r)"''"' 87re25^(r)7' 



j(q) = -87-A2(T)a(q) ^ -^-^^a(q) A^T) = (95) 



The set of equations (|93l94p allows to generalize London equation. From the equation we can find the value of 
the gap: 

^ 0, 2^ (96) 

2dQ + 27 (q2 + i^a2) 

Then the equation for current has the form: 

j(q) = -87-- -3a(q)^Q(q,a)a(q) (97) 



(2do + 27 (q2 



4e2 2 



This equation is the nonlocal and nonlinear generalization of London equation in the long wavelength limit g — > 0, 
because the kernel Q is the function of q and of the magnetic field a(q). However the equation is correct when 
the magnetic field is much weaker than the critical field H <^ He- 

Let's neglect by the nonlinearity, that is we suppose that the kernel Q is function of q only. Then we have: 

- -87^ ,,/°^f ,,2 a(q) = -87^-^^^a(q) « -87^^^^^a(q) ^ Q(g)a(q), (98) 
c (2do + 27q^) (l + d~^j ^ [L + L^q ) 
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where Igq <C 1, ^ = = Iq is the coherent length at temperature T = 0. Thus we have the nonlocal kernel 

Q{q), where radius of a nonlocality is equal to the coherent length Iq. This result corresponds to nonlocal Pippard 
electrodynamics (m the long wave limit). This fact proves a nonlocality of the obtained functional of free energy of a 
superconductor (|79[) . For generalization in case of the large value IqQ it is necessary to expand the free energy (|79p in 
degrees of q. 

V. CONCLUSION. 

In this paper we developed the microscopic approach for a finding of the functional of free energy of a superconductor 
with help the direct calculation of the vacuum amplitude. This value is calculated on the dressed propagators, which 
takes into account the interaction of a free fermion with the fluctuations of pairing (with the condensate of pairs). 
As a result of such interaction the dispersion law of quasi-particles is changed and anomalous propagators appear. 
This means, that a spontaneous symmetry breakdown takes place. After consideration of interaction of particles with 
the fluctuations of pairing all characteristics of a system must be calculated over the new vacuum with the broken 
symmetry. As a result, the free energy is function of the amplitudes fl = ri(AA+), and their observed value minimizes 
the free energy. 

With help of the developed microscopic approach we have obtained the free energy of a spatially inhomogeneous 
superconductor (the system with a spatially inhomogeneous distribution of the order parameter A(r)) for all 
temperatures. The cause of the nonhomogeneity is the dependence of the interaction constant of electrons on 
coordinates and the influence of a magnetic field. The obtained expression for the free energy is valid for arbitrary 
spatial variations of the order parameter. In the long-wave limit the functional is similar to Ginzburg-Landau 
expansion, but the coefficients of expansion depend on coordinates. Moreover, the obtained expression shows the 
nonlocality of the order parameter to the extent that the value of the energy gap in a point is determined by the 
distribution of the energy gap in some neighborhood of this point. 

Above-mentioned results have been generalized on the case of a presence of a magnetic field and a current. The 
analysis of the obtained functional of free energy ([7^ shows, that the nonlocality of a magnetic response of a 
superconductor is the result of the nonlocality of the order parameter. The equations of superconductor's state are 
extremals of the functional obtained by variation over the gap A and the magnetic field a. In the high-temperature 
limit the obtained equations is analogous to Ginzburg-Landau equations. In the low-temperature limit the equations 
show the nonlocal connection between the magnetic field and the current according to the Pippard low. 
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